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Oh. Abstract: Cold electroweak baryogenesis was proposed as a scenario to bypass generic 

r~| ! problems of electroweak baryogenesis within the Standard Model. In this scenario, baryo- 

genesis takes place during an electroweak symmetry breaking transition, which is also 
'k/ \ responsible for preheating after inflation. In the simplest modelling of the scenario, only 

;h ' two parameters remain undetermined: The Higgs mass and the strength of CP violation. 

Using full real-time lattice simulations, we compute the dependence of the asymmetry on 
these parameters. 
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1. Introduction 

Any scenario of baryogenesis aims at reproducing the observed matter-antimatter asym- 
metry of the Universe, usually quantified as the baryon-to-photon number density ratio 

i, 

^ = 6.1x10-10. (1.1) 

In the original electroweak baryogenesis scenario, this is achieved using Standard Model 
(SM) physics (or electroweak physics of a supersymmetric extension of the SM), at a first 
order electroweak phase transition Q (see |y] for a review). Within the SM proper, the 
electroweak phase transition is a cross-over and cannot accommodate the required out- 
of-equilibrium conditions for successful baryogenesis Q. Furthermore, (at least at finite 
temperature) the SM CP-violation in the CKM fermion mass matrix is insufficient by 
many orders of magnitude (5|, ^. In Cold Electroweak Baryogenesis, a period of infiation 
is assumed to end at the electroweak scale j^, ^, and subsequently electroweak symmetry 
breaking takes place at zero temperature, is strongly out of equilibrium through the process 
of tachyonic preheating, and is also responsible for the (re)heating of the Universe. Further 
details on different aspects of the scenario can be found in y, |l^, |l^] (low-scale inflation), 
fl^ , |l^ , |l4| (electroweak tachyonic preheating), |l^, 16, |l^ (generation of the asymmetry). 



|l8f| (SM CP-violation at zero temperature). This work is a continuation of ||17|], where full 
lattice simulations of the electroweak transition including CP-violation were carried out. 



2. The SU(2)-Higgs model with CP-violation 

We study the model described by the action: 

(2.1) 
where e is such that the energy density in the ground state is zero, //eff is a time-dependent 
effective mass for the Higgs field and n parametrises the strength of effective CP violation. 
In Cold Electroweak Baryogenesis, Higgs symmetry breaking is triggered by a coupling to 



an inflaton. In Inverted Hybrid Inflation Q, e.g. the one in [11 1 



filsmU = [/^' - A^<^a2(i)] </>t^, (2.2) 

where a{t) is the time-dependent expectation value of the inflaton field. As in [^ we will 
specialise to the instantaneous quench, 

^,lf,{t < 0) = fi^ A^2ff(i>0) = V, e = ^V(4A). (2.3) 

The case of non-zero quench time will be treated in a separate publication |1S]. A suf- 
ficiently rapid change of sign in //^g induces a spinodal instability with large occupation 



numbers, enabling us to use a classical approximation ||13|, |15|| . 

The CP-violating term is to be thought of as a lowest dimensional effective contribution 
from a theory beyond the SM ||5|, ^. The parameter k is dimensionful and can be written 
in terms of a dimensionless parameter 6cp as 






(2.4) 



(2.5) 



For definiteness we have used mw as the mass scale in (^), mw = 9v/2, with v the 
vacuum expectation of the Higgs field, v'^ = fi'^/X- Experimentally, mw ^ 80 GeV and 
V = 246 GeV, which fixes the gauge coupling to be 5 ~ 0.65. We use g = 2/3 corresponding 
to mw = 82 GeV. We will allow the Higgs mass to vary relative to the W mass, determining 
the Higgs self-coupling through 

m,H \ 8A 

^mwj g 
For the cases considered here, m?^/m,'y^ = 2, 3, 4. 

We choose 5cp to be in the interval (5cp = [0, 1]. The aim is to interpolate to very small 
values of (5cp since, as we will see, to reproduce the observed baryon asymmetry we will 
need 5cp = 0(10^^). Ideally, we are looking for a linear regime at small 5cp- In fl^l) ^^ 
found that the dependence is non-linear for the range of 5cp used there, and the present 
work zooms in on the interval between zero and the first non-zero value of (5cp in |jl^ . 

2.1 Equations of motion and observables 

The action is discretised on the lattice and the classical equations of motion are derived 
(for details on the lattice implementation, see [|l^]). In the continuum, they read: 

'd^D^^ + iils{t) - 2A</.V - kTiF'^^F^,) </> = 0, (2.6) 



Do (^El - 2K(PUBi] - eki^Di (^B''^ + 2k<pUE'^) + Jt = 0- (2.7) 

with E^ = F^Q, Bl = ekimEi'^/2, Di is the adjoint covariant derivative D'jf = S^'^d^+eabc^'li 
and j" is the SU{2) current from the Higgs field, 



^a 



j- = i{D^<P)^^^-i^^^D^^. (2.8) 



The Gauss constraint, 



Dk{i^E-,-2KcpUBi]+JS = 0, 



should be imposed on the initial condition, and will then be conserved by the equations of 
motion. 



As was mentioned in [17|, a good lattice implementation of FF results in implicit 
equations of motion, which require iterative solving. In combination with the number of 
terms arising from FF, this means an increase of computer running time by roughly a 
factor 10. The work presented here amounts to about 30 CPU-years. 

We use periodic boundary conditions with spatial volume L^, and study the evolution 
of the Higgs expectation value. 



1 

24^2 



A^w = :r— ^ / d x ej^feTr 



iTr$t$' ^ 



^~^ = ^yd^x^, (2.10) 

the Chern-Simons number, 

Ncs{t) - A^cs(O) = ^ j dt j d'xTiF^'^F^,, (2.11) 

and the Higgs winding number, 

{d^u)u\djU)u\dkU)u^ 

(2.12) 

A^w is integer and can only change if there is a zero of the Higgs field. Such a zero 
is energetically unfavourable. Once the Higgs field has settled near the bottom of the 
potential (^^ « 1) and the temperature is relatively low (« 50 GeV ||l^), no further 
changes should be seen in Higgs winding. In equilibrium, at sufficiently high temperature, 
winding number changing transitions occur when the Higgs field goes through sphaleron- 
like configurations. (Way) out of equilibrium and in the presence of a large number of Higgs 
zero's, winding number can change readily. Winding and unwinding during the tachyonic 
electroweak transition was studied in [^, |2^ . 

Acs is integer in the gauge vacua, and equal to Aw- At finite temperature (or finite 
energy density, out of equilibrium), it need not be integer and may be very different from 
A^w Still, as the system thermalises to a low temperature, we would expect the Chern- 
Simons number to relax to a value close to the winding number. This is what we see 
happening for late times. We will use this fact to simulate only until the Higgs winding 
has settled and the transition is over. This winding will then tell us what the asymmetry 
in Chern-Simons number would be, had we waited for it to settle. 



2.2 CP symmetric initial conditions 



We initialise our Higgs field using the "just the half" method |T5| (also used in |21, 22, 23 



12, |T^). An ensemble of configurations is generated reproducing the quantum two-point 
functions in the vacuum before the quench. 

^'^'^'^'^^ = 2;^;:' ^^i'^k) = Y' ^^ = Vt^^+k^ (2.13) 

with ^k = L~'^''^ f d^x e~ (p{'K), and similarly for tt. Gauge fields Ai are zero initially, 
with their canonical momenta Ei determined through the Gauss constraint. 

The ensemble of initial configurations is CP-symmetric. However, in a numerical simu- 
lation one only has a finite number of initial configurations available. Let N = j^ ^7=1 ^j^ 
be the numerical estimate for {N) , where N = Ncs or N^ and M is the number of initial 
configurations. Even for 6cp = 0, A^ is typically non-zero because of statistical fluctuations. 
In a plot of N versus 6cp this leads to large uncertainties in the slope dN/d5cp near the 
origin. Previously |17| we dealt with this problem by using the same series of pseudo ran- 



dom numbers for 5cp = and 6cp 7^ 0. Here we avoid it by including the CP-conjugate 
configuration with every randomly generated initial configuration. 
We define the observables 

^JVcs 2 ' ^" ~ 2 ' ^ -^ 

where A^, A^*^^ correspond to the values for a CP-conjugate pair of initial configurations, 
as in ( p. 21 ). Taking ensemble averages, we obviously have 




(A^^J = (ATe,), (A^J = (Ar„). (2.15) 

These observables have the advantage that they cancel out some of the statistical noise. In 
particular. Aw takes integer and half-integer values, which reduces fluctuations. Standard 
errors calculated in terms of A are smaller than for A^. We have 

-i = ( ( tl^^j j - (^.1^^^ = 1 (,^ + l^NNcp) - {N){Ncp)) , (2.16) 

where we used 

a% = {N') - {Nf = {Nip) - {NcP?. (2.17) 

This gives an error estimate of 

2 \{al + {NNcp)-{N){Ncp)) ,_^, 

errors = ^^^^-^ , (2.18) 

where M/2 is the number of pairs of configurations. It may be compared with using simply 
the observable A^ with M random initial configurations, 

2 
evTorl = -^^. (2.19) 



In our case, N ~ —Nqp, in which case the cross correlator in ( ^.18 ) is large and negative 



{{Ncp) = i^)), reducing the error. In the limit of no cross-correlation and M ^ 1, the 
two error estimates ( p.l8| ), ( |2.19 ) coincide. 



Let a prime denote the operation of CP conjugation, 

0'(x, t) = (/.'(-x, ty, A'fc(x, t) = -Afc(-x, tf, (2.20) 

where * denotes complex conjugation and T denotes transposition. If (p, A are a solution 
of the equations of motion with ^cp; then (j)', A' are a solution with — (5cp and CP-conjugate 
initial conditions. As before, let N {6cp) denotes the final N resulting from CP-conjugate 
initial conditions without changing 6cp- Since N is odd under CP, it follows that 

N^'^iScp) = -N{-6,p), (2.21) 



as illustrated by the following diagram {S = 6cp) 

</>(x,0) ^Ni6) 
i CP (2.22) 

Expansion in S, 



(A'(x,0) ^ N'{-6) = -N{-6) 



N = No + Ni5 + N25^ + 0{5^), A = Ni5 + 0{5^), (2.23) 



gives 



(N) = {Ni)d+{N2)S^ + ---, (2.24) 

(^N = (^o) + 2{NoNi)5 + (2iVoA^2 + Nf)6^ + • • • , (2.25) 

al = {Nf)6^ + 0{6^). (2.26) 

Note that the zeroth and first order terms are absent in o"^, suggesting a strong reduction 
in statistical noise for small 6cp- 

Because of ( p.21| ) we just need to run with =b(5cp, rather than the actual CP-conjugate 
configurations. We checked this numerically. In the following we will no longer distinguish 
between the exact (A^) and the numerical estimate N. 

3. Numerical results 

In figure 0(left) we show the evolution from a single initial configuration, evolved with 
Sep = ±1 and niH = 2miy. For the two trajectories the Higgs field c/)^ performs symmetry 
breaking in an identical way, settling near 1. At the same time, Chern-Simons number 
grows in an almost symmetric way. Shown here is A^cs('5cp = 1) and — A^cs('^cp = —1) = 
N^^{Scp = 1). The Higgs winding number is truly symmetric after settling near time 10 
(the glitches of magnitude less than one are discretisation errors). Obviously, this is a 
configuration pair with no generated asymmetry, Ajv^ = 0. Figure ^(right) is a similar 
pair of trajectories, but now A^r^ = 1/2. 




Figure 1: Example of CP conjugate pairs of configurations. Shown is the time evolution of (fy^ 
(black), A^cs (red) and N^ (blue) for (5cp = 1, niH /raw = 2. In one case, the CP violation has very 
little effect (left) and /S.m„ — 0, in the other (right) the result is a net difference of /S.m„ ~ 1- 
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Figure 2: The final values of A^w for 96 pairs of configurations. Black is N^{5c-p): red 
Nvi{—5c-[>) = ^N^^ {5cp)- When they are not on top of each other, a net asymmetry has been 
produced; niH/mw = 2, Sep = 1. 



We run until mnt = 100. Figure ^ shows the resulting values of N^ for +(5cp (black) 
and —6cp (red), using an ensemble of 96 pairs of configurations. {^n„) is (a half times) 
the black minus the red values, averaged over the 96 pairs. 

From now on all results are for the ensemble averaged quantities {N^), (Ncs)- Although 
the final asymmetry is what we are ultimately interested in, the full time evolution shows 
complicated features. There is a linear regime during the first rolling off of the Higgs field. 
Then a non-linear back-reaction regime, where the behaviour of (A'cs) can be described 
approximately in terms of diffusion under a time-dependent chemical potential |^, |l^. A 
nice aspect of this description is that the dependence on the CP-violation is clearly linear. 
This regime ends when (A^w) begins moving away from zero, to settle near the final value. 
The change of {N^) is determined by the generated {Ncs) as well as the availability of zeros 
of the Higgs field, and is as such a very complicated process. 

In the following we will express time in units of the Higgs mass, r = niHt. 



3.1 Initial rise 



In the initial tachyonic instabil- 
ity of the Higgs field, low momentum 



modes dominate. In |17] we solved 
for the early time evolution in the 
linear regime, making a homogene- 
ity approximation and treating the 
CP violation as a perturbation. The 
result was that the generated asym- 
metry during the first roll-off of the 
Higgs field is given by: 



(iVo 



647r4(l-hc)2 m 




0.001 r 



0.0001 



(3.1) 



Figure 3: Comparing the initial bump from the full 
simulation (red) to the estimate from the linearised and 
homogeneous equations of motion, cq. (|3.lD (green). 



The parameter c is to be extracted 
from the growth of the magnetic field 
^2 oc exp(2V2cr), and B'^ and 4? 

are taken from the simulations. For the case miilmy^ = 2, c ~ 0.62, for mu/myi/ = v2, 
c ~ 0.67. Figure ^ shows the result for {N^s) from eq. ( |3.1| ) (green) and the full simulation 
(red) for the case mH/mw = 2. Notice that the scale is logarithmic. The discrepancy is 
20 percent up to time 5. When scaling the simulations for various 5cp linearly with (5cp, the 
red curves fall on top of each other. 

3.2 Early back-reaction: Asymmetric diffusion 

The value of cfP' in figure |l| indicates that the system becomes non-linear after r ~ 5. We 
may get some insight into the early subsequent behaviour by considering the diffusion of 
Chern-Simons number under the influence of a chemical potential p^ , |7|, ^]. Making an 
approximation in which (frcf) in the CP-violating term in the action (^]^) is replaced by 
its spatial average and making a partial integration, exhibits a time-dependent chemical 
potential-like interaction for Chern-Simons number 

„,2 



d X K 



TiF'^'F, 



^1/ 



dt K — (j) 



~^d 



2— Wir'^Na 
dt 



dt//r^a 



with 



cr(t) 



9^ dt 



(3.2) 



(3.3) 



where we also used the definition (l2.4|) of 5cp and m^^ 



jQ^v'^. The effective diffusion rate 



of Chern-Simons number is 



dt 



{{ND - (iVes)') . 



(3.4) 



""^In equilibrium, V is called the sphaleron rate, which describes the widening of the distribution of Chern- 
Simons number. In the present out-of-equilibrium context, it does not have this straightforward physical 
interpretation. 



The generated Chern-Simons number asymmetry is then deduced to be [0, |l^] 

, ^r(i')r(t') 



(iVc.)(t) = / dt' 

'0 



fcff 



(3.5) 



where Tefj can be thought of as some effective temperature of the relevant low momentum 
modes. 

Because the gauge fields become large as the transition proceeds, the effective diffusion 
rate grows in time. Performing the integration directly from the time-dependent, numeri- 
cally determined T{t) and 4>'^{t) (figure §) and eq. ( |3.5| ), one reproduces not only the initial 
rise, but also the subsequent dip, resulting in an asymmetry with the opposite sign from 
the initial rise (figure ^). Indeed, because the diffusion rate is larger towards the end of 
the transition (r ~ 7), it conspires with the Higgs field oscillation (with negative slope, 
effective chemical potential) to qualitatively change the final outcome. In this argu- 



— d<Ncs >/dm^t 

2 ^ 

— d<(|) >/dm t 




-0.2 



Figure 4: The effective diffusion rate r(t), eq. ( |3.4| ) (red), and the time derivative of the Higgs 
field, proportional to is^^^{t), eq. (|3.3D (black), both in units of niH- 



ment, Tefj is an adjustable parameter, in figure ^ taken to be ~ Sttih- This corresponds 
to Tefj = 1.3 TeV, which is quite large. Once chosen, the semi-quantitative agreement for 
different 6cp is convincing. The evolution of (A'cs) is described by eq. ( |3.5D until r ~ 8 — 10. 
The agreement ends around time r = 10, which is also when the Higgs winding begins 
to grow. Apparently, the linear-response treatment for the Chern-Simons number cannot 
account for the dynamics of winding and unwinding. For this, only the full non-linear 
simulations give a correct picture. 



3.3 Intermediate times: Higgs winding creation 

As we have seen, up to r = 10 the statistical treatment of the Chern-Simons number is 
quite successful; There is a net (Ncs), but still a tiny (A^w)- Energetically, the two are 
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Figure 5: Comparing eq. ( |3.5[ ) (black) for (iVcs) to the full simulation (red), raulniw — 2, (5cp = 
1/8, 1/4, 1/2, 3/4, 1. The red lines are curves for all (5cp, simply rescaled to (5cp — 1. The blue 
lines are (N^), also rescaled. Notice that the average winding number does not move until around 
r = 10, the time of the first minimum of {4P'){t). 



strongly favoured to end up near each other at later times. This means that one has to 
adjust to the other^. 

Higgs winding only changes when 
there is a zero of the Higgs field. The 
average Higgs field 0^ continues to 
oscillate some time after the transi- 
tion (figure]^. When it is low, the 
probability of zeros in (j) itself is high. 
The creation and evolution of (near) 




Figure 6: The time history of (0^) (black), (A^cs) (red) 
and {Nvi) (green) for niH/mw = 2, 6cp ~ I. 



zeros was studied in |20|, where it 
was seen that they indeed act as nu- 
clei for winding number change as 
well as sphaleron-like transitions. It 
was also seen, that multiple "genera- 
tions" of (near) zeros are generated, 
corresponding to subsequent minima 
of the Higgs oscillations. First gen- 
eration nuclei are the most numerous, subsequent generations are less populated. 

The existence of such zeros suggests why in the first Higgs oscillation, around r = 12 
the Higgs winding is able to adjust to the Chern-Simons number (figure ^). For late times 
the winding number can no longer change, except through true sphaleron transition, for 
which the time scale at these temperatures is very long compared to the time scale of the 
simulation. We can estimate it through Tsph oc e~^=p'"^, with E^ph the sphaleron energy of 
order 10 TeV or 60 niM. At time 2::; 1 00m 
^A similar situation has been studied in fca, Gi 
size of the winding and Chern-Simons number 'blobs' is an indicator whether A^w adjusts to A', 



, a Bose-Einstein fit to the particle distribution 
, where it was seen that in single trajectories the relative 



functions gives T/mn — 0.4 |12], suggesting that the sphaleron rate is indeed very smaU at 
these times. It also suggests that Tgff as extracted from the asymmetric diffusion (section 
|3.2| ) should be interpreted with care. At longer times, the Chern-Simons number will settle 
close to the winding number value^. 

3.4 Dependence on Higgs mass 

The end result turns out to be very sensitive to the Higgs to W mass ratio. Here we present 
results for niH = v2?n.vK and \/?>mw- In the former case (figure 0, left), the overall sign 
is opposite to what we saw in figure |6[ In the latter (figure j^, right), we are apparently in 
an intermediate case, where although there is still the initial linear regime, the dynamics 
conspires to give a final asymmetry consistent with zero. The equation of motion of the 
Higgs field depends on the time derivative of the Chern-Simons number, and the frequencies 
and phases of these oscillation can conspire to give asymmetries of opposite signs. 
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Figure 7: The time history of ((/>^) (black), (iVcs) (red) and (TVw) (green) for different mass ratios; 
mH/mw — V2 (left) and V3 (right); <5cp = 1. 

The large dependence on the Higgs mass is reminiscent of the situation for the anal- 
ogous Abelian-Higgs system in 1-1-1 dimensions |jl^. There, we were able to span a much 
larger range of masses, and the resulting curve looked quite complicated (figure ^ (left)). 
For details about similarities and differences between the two studies, see [|l5[. At present 
we do not have the computational resources to perform an equally thorough study in 3-1-1 
dimensions. But using the result for rriH = rriy/ from ||T^ to guide us^, we can tentatively 
draw a plot of the mass dependence, figure ^ (right). 

3.5 Dependence on CP-violation 

The nonlinear behaviour at intermediate times might also destroy the linear dependence of 
the final asymmetry on (5cp. To study this we vary 5cp using the same initial configurations 
for all (5cp. To get meaningful errors in the case ran /raw = 2, we had to increase the number 



^This we have checked for a few configurations, running to r = 500. 

*An error in the apphcation of 5cp in [|17|] has been corrected. Furthermore, results of |l7|] should be 
multiplied by g^ = 4/9, the initial conditions for the case mn = rnw at small k were different ('thermal', 
resulting in somewhat smaller final results), and the quoted value refers to (A^cs) rather than {N„). 
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Figure 8: The Higgs mass dependence of the asymmetry in the analogous model in 1+1 (left) and 
3+1 (right) dimensions; ^cp = 1, k is the analogue of Jcp- Lefthand plot from ||l^. In the righthand 
plot, blue points are the simulation presented here for Sep — 1, red are the results from the fits to 



the (5cp-dependence (see below) and black is the mH/inw = 1 result from |17 
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Figure 9: Final (A^w) as a function of Jcp. Left: mH/mw — 2, right: iuh /niw — "\/2- The full 
red line is a linear fit, the dotted lines represent ±lcr in the fitted slope. 

of CP-conjugate pairs to 192. Figure ^ shows the final value of the average winding number 
vs. Sep up to 5cp = 1- Within errors, the dependence is consistent with linear. The fits in 
figure ^ and the one final value from figure (right) lead to an asymmetry 

(iV„) = (0.075 ± 0.006)(5cp, mn = V2mw, 

{N„) = (0.005 ± 0.020)5cp, mH = V3mw, 

(N^) = (-0.0359 ± 0.0040)(5cp, mn = 2mw. (3.6) 



4. Conclusion 

Given the final ensemble average of the winding number, we can make an estimate for the 
generated baryon asymmetry. We use 
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The entropy s is given in terms of the reheating temperature T^eh and g^: the number 
of relativistic degrees of freedom and the reheating temperature deduced from the initial 
energy density in the Higgs potential, T/mn — 0.45. We also assume, as discussed earlier, 
that the late times (A'^cs) will be equal to (A^w) at the end of our simulation. We have 

(4.2) 
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(4.3) 



Compared to [17| the CP-symmetric initial conditions and the emphasis on the Higgs 
winding number allowed us to get a much clearer signal without much larger statistics. 
This was necessary in order to zoom in on the range of 5cp where the dependence is linear. 
In particular, we were able to pin-point the time at which the asymmetry is generated to 
the first minimum of the Higgs field evolution. This is when the average winding number 
is able to change and accommodate the initial asymmetry in the Chern-Simons number. 
In the range 5cp = [0, 1], the asymmetry is linear in 5cp, allowing us to interpolate to the 
very small values relevant for the observed asymmetry. To reproduce the observations 
( |1.1| ), we need (5cp — 2 x 10^^ (m// = 2mw)- Presumably, 6cp should be somewhat larger 
than this, when taking into account the dynamics of the inflaton, fermions and additional 
gauge fields, which may in various ways affect the dynamics of the SC/(2)-Higgs system. 
In particular, the assumption of an instantaneous quench leads to quite wild behaviour. 
In very slow quenches, the system may never be sufficiently out of equilibrium, and the 
asymmetry should be correspondingly small. The dependence on the quench time will be 
presented in a separate publication |jT^ . 

The mass of the Higgs field in the Standard Model is expected to be smaller than 200 
GeV ~ 2.5 mw- We have probed the allowed region and found a dramatic dependence on 
niH- Whether or not this effect survives at finite quench times is not yet known, and it 
is clear from the semi-analytic linear treatment in section ^^, that the generic sign of the 
asymmetry is opposite to that of Sep for mn = 2 m\Y . Still, both for mn = V^ mw and 
for rriH = rnv/ |l^ the final result has the opposite sign, i.e. the same sign as 5cp- 
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